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Abstract 

We have investigated the role of locahzed d bands in the dynamical response 
of Au, on the basis of ah initio pseudopotential calculations. The density- 
response function has been evaluated in the random-phase approximation. For 
small momenta, we have found a double peak structure in the energy-loss func- 
tion, which results from the presence of d electrons. These results are in good 
agreement with the experimentally determined optical response of gold. We also 
analyze the dependence of the dynamical structure factor on the wave vector q. 



1 Introduction 

Many experimental and theoretical work have focused on noble metal systems 

in order to get a better understanding of the role played by localized d electrons on 
the properties of delocalized, free-electron-like, electrons. 

The behavior of noble metals has been considered by Ehrenreich and Philipp 
These authors used the data obtained in experimental measurements of the reflectivity 
from Cu and Ag @] to obtain the optical spectra. They separated the Drude and 
interband contributions to the optical response, explaining why the energy of the 
plasma resonance in Ag is shifted down in energy. In a later work. Cooper, Ehrenreich 
and Philipp ||^ extended this analysis to Au. 

Gold and copper present no decoupling between sp and d orbitals, and a com- 
bined description of these one-electron states is needed to address both structural and 
electronic properties of these materials. 

Due to the rapid progress in computing technology, realistic calculations of the 
density-response function of solids have become feasible in the past few years. Re- 
cently, the dielectric matrix as well as the energy-loss spectrum have been evaluated 
numerically within a first principles framework for several simple (Al |Q, Cs [0], Be [H 
) and noble metals, as well as for some insulators like Si 0, LiF [T^ or AI2O3 



In the case of noble metals, the energy-loss function for Cu has been evaluated 
in Ref. 4 and plasmonic excitations in Ag have been reported in Ref. 5. Within the 
same many-body framework, ab initio calculations of the electronic stopping power of 
real solids |jl2[ and excited electron lifetimes |]13| have also been reported. 



In this paper we report an ab initio evaluation of the dynamical response function 
of Au as computed in the random-phase approximation (RPA) after an expansion of 
all and 5(i^° one-electron states in a plane wave basis. 



2 Formalism 

The linear density-response function x(r, r'; u) of an interacting electronic system is 
defined by the equation 

r\r, u) = J dr' x(r, r'; V'^\r', u), (1) 
where p*"'^(r, cu) is the electron density induced by an external potential V^*^^*(r', cu) 



Tl 



For periodic crystals we Fourier transform x(r, r';c(j) into a matrix Xg,g'(^!^) 
which, in the framework of time-dependent density-functional theory (TDDFT) [|T5[| , 
satisfies the matrix equation 



XG,G'(q, ^) = XG,G'(q. ^) + EE XG,G"(q, ^) 

G" G'" 

^^G"(q)'^G",G"' + ^G",G"'(q'^)] XG"',G'(q'^)- (2) 



The wave vector q is in the first Brillouin zone (BZ), G and G' are reciprocal lattice 
vectors, fG(q) = 47r/|q-|-Gp are the Fourier coefficients of the bare Coulomb potential, 
the kernel i^Q Q,(q, cj) accounts for short-range exchange-correlation effects (in the 
RPA, i^'QQ,(q, u;) = 0) and Xgg'I^'^) the Fourier coefficients of the density- 
response function of noninteracting Kohn-Sham electrons: 

/ \ 1 /k.n /k+q,n' 

XG,G'(q>^) = o EE 



^ k ^k,n - ^k+q,n' + + ^V) 

X(0k,n|e^'(''+^)-l0k+q,n')(0k+q,n'|e'(''+^')-'-|0k,n), (3) 

where the second sum runs over the band structure for each wave vector k in the first 
BZ, /k,n represents the Fermi-Dirac distribution function, rj is a. positive infinitesimal, 
and fl represents the normalization volume. 0k,ra(r, cj) and i^k.n are the Bloch eigen- 
functions and eigenvalues of the Kohn-Sham Hamiltonian of density-functional theory 
(DFT) 

The inverse dielectric function is connected to the density-response function by the 
following relation: 

eG!G'(q>^) = (^G.G' + ^^G'(q)XG,G'(q,tu). (4) 

The properties of the long-wavelength limit (q 0) of the dynamical density-response 
function are accessible by measurements of the optical absorption, through the imaginary 
part of the dielectric response function eG=o,G'=o(q = O,^^). 

For the evaluation of the one-electron Bloch states, we first expand them in a 
plane-wave basis, 

</>k,n(r) = 47E^k,n(G)e^(''+«)-, (5) 

V V G 



with kinetic-energy cutoff of 75 Ry (~ 10^ G- vectors) and then evaluate the Kohn- 
Sham equation of DFT with a full description of the electron-ion interaction based on 
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Figure 1: Convergence for the energy-loss function for different number of bands and 
samplings for q = 0.16(1,0,0) x (27r)/a (a = 3.01 A is the experimental lattice constant) 
and G = 0. Optical data are from Ref. 19. "b" stands for bands and "x" for the sampling 
(See the text for reference). 

the use of a non-local, norm-conserving ionic pseudopotential Finally, we evaluate 
from Eq. d^) the Fourier coefficients Xcci^y^): solve the matrix equation, Eq. 
(^, for the Fourier coefficients of the interacting density-response function, which we 
evaluate in the RPA. 

We have performed calculations for the energy-loss function using 25, 30 and 40 
bands, and samplings over the BZ on 10 x 10 x 10, 12 x 12 x 12 and 14 x 14 x 14 
Monkhorst-Pack meshes [0. 25 bands and a 12 x 12 x 12 sampling turned out to be 
enough to reach convergence (See Fig. |l|), but 30 bands were used to guarantee good 
values for the highest energies under study. 



3 Results 

Fig. shows our ab initio calculations (solid line) for both real and imaginary parts 
of the eG=o,G'=o(ct = 0,uj) dielectric function of gold for a small momentum transfer 
of |q-|- G| = 0.13 Oq ^ (flo is the Bohr radius) together with the optical data (q = 0) of 
Ref. 19 (dashed line). 

In Fig. ^the corresponding values of the so-called energy- loss function Im — Cg^gI^j ^) 
are displayed. Our results, obtained for a small but finite momentum transfer, are in 
excellent agreement with the measured dielectric function, both showing a double peak 
structure in the energy-loss function. The solid line represents the solution when we 




Figure 2: Real and imaginary parts of the eG,G(q)'^) dielectric function of Au, for q = 
0.16(1, 0, 0) X (27r) I a (a = 3.01 A is the experimental lattice constant) and G = 0. Solid and 
dashed lines represent our calculations and the optical data of Ref. 19, respectively. 



consider the crystalline local-field effects appearing through the dependence of the 
diagonal elements of the interacting response matrix xq,^q,i(o^^iS) on the off-diagonal 
elements of the polarizability XQQ,(q, cj). The dashed-dotted line describes the be- 
havior when this effect is not taken into account. Even if the shape of both curves 
is quite similar, the solid line accounts better for the experimental result than the 
dashed-dotted line, the difference being bigger for higher energies. The small finite 
q- value used keeps all the features of the optical spectrum and therefore can be directly 
compared with the measured one. 

In order to investigate the role of localized d bands in the dynamical response of 
gold. Fig. ^ also shows (dotted line) the energy-loss function calculated from an ah 
initio pseudopotential with the 5ci shell assigned to the core. The plot shows that 
a combined description of both localized hd bands and delocalized electrons is 
needed to correctly address the actual electronic response of gold. 

For Au, as well as for the other noble metals, one expects to observe a resonance 
near 9 eV corresponding to the density of free electrons. The role of the Au d bands is 
to provide a polarizable background which lowers this free electron plasma frequency. 
We note from Fig. that near 5 eV the real part of the dielectric function (Re e) is 
zero; however, the imaginary part (Ime) is not small, due to the existence of interband 
transitions at these energies which damp the free-electron plasmon. 

d-bands are also responsible, at higher energies, for a double peak structure in the 
energy-loss function, which arises from a combination of band-structure effects and the 
building up of (i-like collective modes. Since these peaks occur at energies (~ 25 eV 
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Figure 3: As in Fig. ^ for the energy- loss function of Au. The dotted line represents the 
result of assigning the 5d shell to the core and the dashed-dotted line the result of excluding 
the crystalline local- field effects. 

and ~ 35 eV) where Re e is nearly zero (see Fig. |^) they are in the nature of collective 
excitations, the small but finite value of Ime at these energies accounting for the width 
of the peaks. A similar double-peak structure has been obtained for the other noble 
metals and seems to be a general feature stemming from the particular (i-electron band 
structure (localization and hybridization with the delocalized s-electrons) of the three 
noble metals. 

In order to see the dependence of the energy-loss function on the momentum 
transfer q + G, we show in Fig. ^ the RPA dynamic structure factor for G = 0, 
as obtained for various values of q along the (100) direction. As long as the 5d shell 
of Au is assigned to the core, one finds a well-defined free-electron plasmon for wave 
vectors up to the critical momentum transfer where the plasmon excitation enters the 
continuum of intraband particle-hole excitations. This free-electron plasmon, which 
is known to show a characteristic positive dispersion with wave vector, is found to be 
completely damped when a realistic description of 5d orbitals is included in the calcu- 
lations. At higher energies and small momenta (q < 0.7), rf-like collective excitations 
originate a double-peak structure which presents no dispersion. For larger values of 
the momentum transfer q + G, single-particle excitations take over the collective ones 
up to the point that above a given cutoff, the spectra is completely dominated by the 
kinetic-energy term as expected. 
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Figure 4: RPA energy-loss function of Au along the (100) direction, for various values of 
the momentum transfer |q-|- G|: 0.16, 0.33, 0.66 and 1.33, in units of 27r/a (a = 3.0lA). 

4 Conclusions 

We have presented ah initio pseudopotential calculations of the dynamical density- 
response function of Au, by including rf-electrons as part of the valence complex. In 
the long- wavelength limit (q — 0), rf-bands provide a polarizable background that 
lowers the free-electron plasma frequency, rf-electrons are also responsible for a full 
damping of this s-like collective excitation and for the appearance of a d-like double- 
peak structure in the energy-loss function, in agreement with the experimentally de- 
termined optical response of gold. We have analyzed the dependence of the dynamical 
structure factor on the momentum transfer, and we have found that, for values of the 
momentum transfer over a given cutoff, a less-pronounced double-hump is originated 
by the existence of interband electron-hole excitations. The dynamical response func- 
tion in Cu and Au shows very similar features, although in Au the peaks appear for 
higher energies. 

We acknowledge partial support by the Basque Hezkuntza, Unibertsitate eta Iker- 
keta Saila. 
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